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Abstract 

Shell effects in the coordinate space can be seen with degenerate Fermi vapors 
in non-uniform trapping potentials. In particular, below the Fermi temper- 
ature, the density profile of a Fermi gas in a confining harmonic potential 
is characterized by several local maxima. This effect is enhanced for "magic 
numbers" of particles and in quasi-lD (cigar-shaped) configurations. In the 
case of a multi-component Fermi vapor, the separation of Fermi components 
in different spatial shells (phase-separation) depends on temperature, number 
of particles and scattering length. We derive analytical formulas, based on bi- 
furcation theory, for the critical density of Fermions and the critical chemical 
potential, which give rise to the phase-separation. 
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I. INTRODUCTION 



The Fermi quantum degeneracy with trapped dilute vapors of 40 K atoms has been experi- 
mentally achieved in 1999. [1J The s-wave scattering between Fermions in the same hyperfine 
state is inhibited due the Pauli principle. Moreover, for 40 K atoms the p-wave scattering 
length becomes negligible in the limit of zero temperature. It follows that at very-low tem- 
peratures the dilute Fermi gas, in a fixed hyperfine state, is practically ideal. Nevertheless, 
the effect of interaction could be very effective for a Fermi vapor with two or more hyperfine 
states. In this paper we consider a Fermi vapor confined in a harmonic potential that models 
the trap of recent experiment. Q These trapped Fermi gases are quite interesting because 
the quantum degeneracy shows up not only in momentum space, as in uniform systems, but 
also in coordinate space. 



II. IDEAL FERMI GAS AT FINITE TEMPERATURE 

Let us consider an ideal Fermi gas in external potential described by the non-relativistic 
field 

^( r ) = 5l0 a (r)a Q , (1) 

a 

where Q (r) = (r|a) is the single-particle eigenfunction with eigenvalue e a and a a is the 
lowering Fermi operator of the single-particle eigenstate |a). The grand canonical thermal 
average of the Fermi spatial density ip + (r)^)(v) is given by 

r»(r)H(^)>) = E jtS + i ' (2) 

where fi is the chemical potential and (3 = 1/ (kT) with k the Boltzmann constant and T the 
absolute temperature. @ Note that the thermal average of a generic operator A is defined 



as 



Tr[Ae~^~^] 



where H = J2 a e a( 1 t^ 1 a is the Hamiltonian of the system and N = J2 a ®i is the number 
operator. The average number N of particles of the system, given by 

a a ° ~ x 

fixes the chemical potential. 

In the case of a harmonic external potential U(r) = (m/2)(uj\x 2 + u^y 2 + ^Iz 2 )-, one 
finds the Fermi density profile by using the Eq. (2) and the eigenfunctions 4> ni n2n :i 

(r) of the 

harmonic oscillator 



V ; ^ e j3R(w 1 (ni+l/2)+W2(n2+l/2)+W3(n 3 +l/2)-At) _|_ 1 " V ; 

niri2".3=0 



Because the Fermi gas is ideal, one has ni „ 2n3 (r) = (f> ni (x)(f) n2 (y)(f> n3 (z), where </> n (x) is the 
eigenfuction of a ID harmonic oscillator with frequency u and quantum number n. This 
eigenfunction can be found by means of the recursion relation || 

<p n (x) = —= \f2ax<p n -i(x) - \Jn- 10 n _ 2 (V) , (6) 

where 0o( a; ) — a" 1 / 2 ?!- 1 / 4 ^ " 2 ^ 2 / 2 and 4>\{x) = -\/2ax(f) (x) , with a = (mw jK) 1 ! 2 . 
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FIGURES 




X 



FIG. 1. Density profiles for an ideal Fermi gas in a ID harmonic trap. Temperature: 
kT/fiu = 10 -3 ; at this temperature the results coincide with zero-temperature ones. ./V is the 
number of Fermions. Lengths in units an = (Ti/mu:) 1 ! 2 and densities in units a H 3 . 

The Eq. (5) and (6) are used to numerically calculate the exact density profile. Such 
exact density profile can be compared with the semiclassical one, that is given by 0] 

n(r) = ^/ 3/2 (e^M>) , (7) 
where A = {2nh 2 (3/m) l l 2 is the thermal length and 

j poo X n ~^ 

/»(*) = =pr / dx , (8) 

T(n) Jo z 1 e x + 1 

with T(x) the factorial function. 
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FIG. 2. Density profiles for an ideal Fermi gas in a ID harmonic trap at finite temperature, 
with N = 20 and 50. N is the number of Fermions. a) kT/huj = 1/10; b) kT/Tiuu = 1/2; c) 
kT/hu! = 1; d) kT/hu) = 5. Units as in Figure 1. 

In the limit of zero temperature, with \i = Ep the Fermi energy, the semi- 
classical spatial distribution is given by the Thomas-Fermi approximation n(r) = 
(2m) 3 / 2 /(6vr 2 ft 3 ) (E F - U(r)) 3/2 (E F - U(r)), where 6 is the Heaviside step function. The 
Eq. (7) is the generalization of well-known formula for an ideal Fermi gas in a box that is ex- 
act in the thermodynamic limit. |2| Note that for \z\ < 1 one has f n (z) = Y^i{~^-) 1+1 z l / i n ■ 
Moreover, by using g n (z) = —f n {—z) instead of f n (z), one finds the spatial distribution of the 
ideal Bose gas in external potential. |||6| The Fermi energy Ep and the Fermi temperature 
Tp are easily obtained by imposing the normalization condition to the Thomas-Fermi distri- 
bution. For a D-dimensional harmonic oscillator with a geometric average u = (ui... ujp)) l l D 
of the D frequencies, one finds the semiclassical result 

E F = kT F = {D\N) 1/D hu . (9) 

For D = 3 one recovers the familiar result Ep = ^^N^-^huj. 

In a recent paper we have performed a comparison between exact and semiclassical 
results. [|j] We have found that when (kT/Tiuj) > 1, where uj = (u \uj 20J -z) 1 ^ ■, there are no 
appreciable deviations between exact and semiclassical results. Instead, for an isotropic 
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harmonic potential, when (kT/hu) < 1 some differences are observable, in particular for 
"magic" numbers of particles (N = (p+l)(p+2)(p+3)/6, where p is a natural number) that 
correspond to a complete shell occupation of single-particle energy levels. The differences 
are reduced by increasing N showing that semiclassical approximation provides an excellent 
representation of Fermi distribution for a wide range of parameters. In correspondence of 
the "magic numbers", the exact spatial density profile shows local maxima, which suggest 
a spatial shell structure. The magic numbers are particularly stable; in fact, for small 
variations of the chemical potential fi the magic number N remains unchanged. The shell 
structure in the density profile is washed out by increasing the number of particles and is 
completely absent in the semiclassical approximation. 

In the quasi-lD case, namely a cigar-shaped gas where uj\ = co> 2 >> &3, the shell effects 
are strongly enhanced. This can be shown by investigating the one-dimensional problem, as 
recently done by Kolomeiski et al || at T — 0. We extend the calculations of Kolomeiski at 
finite temperature. In Fig. 1 we plot the density profile of a ID ideal Fermi gas in harmonic 
potential as a function of the number N of particles at kT/hu = 10 -3 . We have verified that 
at this temperature the density profiles coincide with the zero-temperature ones. The results 
are obtained by numerically evaluating expressions (5) and (6) in the ID case. The local 
maxima, whose number grows with N, are clearly visible for a small number of particles. In 
Fig. 2 we show the density profiles with 20 and 50 particles as a function of temperature. 
Remarkably, the local peaks are no more distinguishable for temperatures well below the 
Fermi temperature Tp = Nhu/k. Thus, to see spatial shell effect on the ID density profile, 
the system should be at temperatures lower than Tp by one or two orders of magnitude. 

III. PHASE SEPARATION WITH TWO FERMI COMPONENTS 

In the recent experiment with dilute 40 K atoms, [|TJ to favor the evaporative cooling, a 40 K 
Fermi vapor in two hyperfine states (|9/2, 9/2) and 1 9/2, 7/2)) has been used. The problem 
of a dilute Fermi vapor with two hyperfine states (components) can be studied by using the 
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s-wave scattering approximation, the mean-field approach and semiclassical formulas. The 
spatial density profiles of the two components of a Fermi vapor can be written as 

n 1 (r) = l/ 3/2 (e^- l/(r) -^ (r)) ) , 
n 2 (r) = l/ s/2 (e^ 2 -^)-Wr))) , (10) 

where ^ is the chemical potential of the i-th component, and g = A-n^a^/ra, with a±2 
the s-wave scattering length between first and second component (a n = a 2 2 = 0). Thus, 
the effect of the second Fermi component on the first component is the appearance of a 
mean-field effective potential. 

When two components have the same number of particles, the onset of phase-separation 
is also an example of spontaneous symmetry breaking. In fact, if the chemical potentials /i« of 
the two components are equal, Eq. (10) always admits a symmetric solution n(r) = rii(r) = 
ri2(r). However, for particle number N larger than a threshold N c the solution bifurcates 
and a pair of symmetry breaking solutions appears. Just beyond threshold the asymptotic 
solutions begin to differ from the symmetric one in a neighborhood of the origin r = 0, i.e. 
at the point of higher density. An analytic formula for the critical chemical potential /i c can 
be obtained by using the bifurcation theory. The system (10) can be re-written at the origin 
as 

F"(x) = , (11) 

where 

x= (x 1 ,x 2 ) = (ni(0),n 2 (0)) 

F"(x) = (iT(x),if(x)) = ( Xl - l/ 3/2 ( e ^-^))),x 2 - l/^e^™)))) . 

For ji < ji c there exist only one symmetric solution of the previous equation, given by 
x* = (x*,x*). At pi — /i c there are two solutions and the Jacobian of the function F M (x) 
must have the determinant equal to zero, i.e. 
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det 



ax 



OlMe^-^) 1 



. 



(12) 



By imposing this condition on Eq. (10) and (11) and setting x* = n(0), one finds 
At T = 0, by the use of the first term of the large z expansion of f 3/2(2), given by 



(13) 



/3/2O 



3^/2 



(\nzf 2 + — (\n Z y 1/2 + ... 
8 



(14) 



one finds analytical expressions for the critical density n c (0) and the critical chemical po- 
tential fi c : 



n c (0) 



7T 



48a 3 



/'c 



24 V a 7 



(15) 



These remarkably simple formulas can be very useful to determine the onset of phase- 
separation in future experiments. Moreover, by knowing the critical chemical potential 
one numerically finds the number of particles. The conclusion is that by increasing the 
interaction between the two components one can use lower number of particles to obtain 
the phase-separation. At finite temperature one can use the first two terms of the large z 
expansion of feniz), i.e. Eq. (14). The analytical formulas at zero and finite temperature 
are in good agreement with the numerical calculations. H 



IV. PHASE SEPARATION WITH MANY FERMI COMPONENTS 

Phase-separation also appears in a Fermi vapor with three or more components. In Ref. 
4 we have numerically solved the three-component extension of Eq. (10) with ayi = a 13 = 
a23 = a. Also for three components the spontaneous symmetry breaking and the phase- 
separation are controlled by scattering length, temperature and total number of particles. 
In particular, one finds that by increasing the scattering length at first one of the components 
separates from the others, which remain still mixed. The separation begins at the center of 
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the trap and, as the scattering length is further increased, also the other two components 
separate (this second phase-separation begins at the interface with the previously separated 
component) and one eventually sees complete phase-separation and the formation of 4 or 5 
shells. 

The Eqs. (15) can be extended to a M-component Fermi vapor with the same number of 
particles in each component. The critical density n c (0) does not depend on the number M 
of Fermi components and one gets the same result of Eq. (15). Instead, the critical chemical 
potential reads fx c = (2M + l)n 2 hu(a H /a) 2 /24. 

V. CONCLUSIONS 

We have analyzed a degenerate Fermi gas in a harmonic external potential. Below the 
Fermi temperature Tp and for a small number of particles, shell effects, like local maxima in 
the spatial density profile for " magic numbers" , are clearly visible. Such effects are strongly 
enhanced in the ID case, but they are visible only for temperatures which are at least one 
order lower than Tp, i.e. far from the semiclassical limit. 

We have also considered a Fermi vapor with many hyperfine states. By using the bi- 
furcation theory, we have shown that the onset of phase-separation appears by increasing 
the scattering length or, for a fixed scattering length, by increasing the number of particles. 
By raising the temperature, a larger scattering length or a larger number of particles is 
needed to obtain the phase-separation. A Fermi vapor with three or more components has 
the same behavior but at first only one of the components separates from the others, which 
remain still mixed. The critical density of Fermions at the origin, which gives rise to the 
phase-separation, satisfies the equation n c (0) = 7r/(48a 3 ), where a is the s-wave scattering 
length. 
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